Universal statistics of non-linear energy transfer in turbulent models 
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A class of shell models for turbulent energy transfer at 
varying the inter-shell separation, A, is investigated. Inter- 
mittent corrections in the continuous limit of infinitely close 
shells (A — > 1) have been measured. Although the model be- 
comes, in this limit, non-intermittent, we found universal as- 
pects of the velocity statistics which can be interpreted in the 
framework of log-poisson distributions, as proposed by She 
and Waymire (1995, Phys. Rev. Lett. 74, 262). We suggest 
that non-universal aspects of intermittency can be adsorbed 
in the parameters describing statistics and properties of the 
most singular structure. On the other hand, universal aspects 
can be found by looking at corrections to the monofractal scal- 
ing of the most singular structure. Connections with similar 
results reported in other shell models investigations and in 
real turbulent flows are discussed. 



One of the basic questions in understanding the physics 
of fully developed three dimensional turbulence is the 
dynamical mechanism characterizing the energy trans- 
fer from large to small scales. According to the Kol- 
mogorov theory (K41) of fully developed turbulence, the 
energy should be transferred downwards from large scale 
to small scales following a self-similar and homogeneous 
process entirely dependent on the energy transfer rate, 
e, and the scale, I. By assuming local homogeneity and 
isotropy, the K41 theory allows us to predict the scaling 
properties of the structure functions S p (l) = ((5v(l)) p ), 
where Sv(l) ~ \v(x + I) — v(x)\. It turns out that 
Sp(l) ~ (e) p/3 l p/3 . The K41 theory has been questioned 
by several authors because of strong scale dependent fluc- 
tuations of the energy dissipation (intermittency). Be- 
cause of intermittency, the scaling properties of the ve- 
locity structure functions acquire anomalous scaling, i.e., 
S p (l) ~ l^ p ', where the scaling exponents ((p) are non 
linear functions of p . 

The universal character of energy transfer statistics has 
been also questioned. Different exponents have been 
measured in anisotropic flows (shear flows || and bound- 
ary layers (J]) and strong dependencies on active quan- 
tities as the temperature in Rayleigh-Benard @j or the 
magnetic field in MHD |@] have been quoted. 
The question of finding unifying and universal aspects 
common to all turbulent flows naturally arises. 
Recently, cascade descriptions based on random multi- 
plicative processes and infinitely divisible distributions 



for random multipliers have been applied in order to clar- 
ify the energy transfer physics. In particular it has been 
recently shown ( pfj[) that a log-Poisson distribution 
is able to provide an extremely good fit to experimental 
data. In JTO^ j authors proposed that the random multipli- 
ers Wl 2 ,L! connecting velocity fluctuations at different 
scales Li and L\ 



6v(L 2 ) = W L2 , Ll Sv(Li) 



(1) 



should follow a log-poisson statistics. 
Log-poissonity should naturally arise in turbulent flows 
as the limit of a Bernoulli fragmentation process between 
infinitely close scales. In authors argued that the 
energy transfer between two adjacent scales l\ and I2 
with log(ix) — log(/2) << 1 can be described in terms 
of only two bricks. The first one corresponds to the 
most singular structure characterized by a local scal- 
ing exponents ho'- Svfa) = (h/h) h °Sv(li). The sec- 
ond brick is a defect-like energy transfer which mod- 
ulates the most singular events by a factor (3 < 1; 
in this second case we would have a typical scaling: 
Svfa) = (3{h/h) h °5v(l{). Let us also assume that, in the 
limit log(/i//2) — * 0, defects along the cascade happen 
with a probability that goes to zero proportionally to the 
logarithm of the scale separation ~ do logilx/lz), where 
do is the parameter which controls how probable a de- 
fect is. It is now simple to show, following [jlo|], that the 
finite-scale-separation transfer must have a log-poisson 
statistics: < {Wl 2 ,l 1 ) p >= (L 2 /L 1 ) h °P- d ^ 1 -i 3 "\ which 
corresponds to the She-Leveque proposal for intermittent 
exponents 0: 



a P ) = h p-d Q (i-pp). 



(2) 



The multifractal interpretation of (@) is that ho is the 
most singular scaling exponent and do corresponds to 
the codimension of the fractal set where the most singu- 
lar scaling is observed. One may argue that the struc- 
ture and the statistics of the most singular event could 
be strongly non-universal. 

As a consequence, a possible scenario can be proposed 
where ho and do could be system-dependent, while (3 
maybe constant in a wider universality class p0| . 
Some evidences of this universal character of (3 have al- 
ready been reported in JTT[ | where it has been shown 
that the differences of intermittent exponents measured 
in Rayleigh-Benard, MHD, shear flows and in boundary 
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layers can all be re-adsorbed by properly changing tig and 
do at constant (3. In |fl| it has been shown that also vis- 
cous effects can be included in a suitable dependency on 
the scale of ho and do , showing for the first time that non- 
trivial intermittent corrections, due to the (3 dependent 
part of ((p) curve, can be detected at viscous scales. 

In the following we will show how log-poisson descrip- 
tion of intermittency and its underlying interpretation 
can be applied for describing some important aspects of 
energy transfer in a class of dynamical models of turbu- 
lence |I^,[l4|,[l3| . In particular, we will be able to ex- 
plain the continuous transition towards the K41 non- 
intermittent statistics in terms of the statistical prop- 
erties of the most singular structure of the model. This 
trend towards K41 statistics is observed in a class of shell 
models at varying (diminishing) the characteristic shell 
separation, i.e. in the so-called continuous limit of in- 
finitely close shells. 

Universal aspects of the shell statistics are recovered by 
measuring the /3-dependent part of the probability distri- 
bution. As proposed by She and Waymire, we find that 
the /3-dependent part of the intermittent statistics is re- 
markably constant for all values of the shell separation A 
explored. 

Shell models have demonstrated to be very useful for 
the understanding of many properties connected to the 
non-linear turbulent energy transfer ( |L3]]- pp|]). The 
most popular shell model, the Gledzer-Ohkitani-Yamada 
(GOY) model ( [^3|- ]2(|), has been shown to predict 
scaling properties for £(p) similar to what is found exper- 
imentally (for a suitable choice of the free parameters). 

The GOY model can be seen as a severe truncation of 
the Navier-Stokes equations: it retains only one complex 
mode representative of all Fourier modes in the 

shell of wave numbers k between fc„ = ko\ n and fc n +i, A 
being an arbitrary scale parameter (A > 1), usually taken 
equal to 2. 

In two recent works [^3|JlJ] the GOY model has been 
generalized in terms of shell variables, it+ and it", trans- 
porting positive and negative helicity, respectively. These 
models have at least one inviscid invariant non-positive 
defined which is very similar to the 3d Navier-Stokes he- 
licity. In the following we will focus on the intermittent 
properties of one of such models at varying the separation 
between shells, A. The time evolution for positive-helicity 
shells is III: 



^ U n — ikn( U n+2 U n+l + ^ U n+l U n-l + CU n-l U n-2) 
— fk^Un + 5 n ,n f + > 



(3) 



and the same, but with all helical signs reversed, holds 
for u~. 

The coefficients b, c are determined imposing inviscid 
conservation of energy, E — J2 n (\u^\ 2 + \u~\ 2 )), and he- 
licity, H = J2 n | 2 ~~ \ u n \ 2 ) ■ I n particular we have: 



b = -(1 + A 2 )/(A 3 + A 2 ) and c = (1 - A)/(A63 + A64). 
Structure functions for these models are naturally defined 

as S p (n) =< (tJ\u+\ 2 + \un\ 2 ) p >~ fc,7 C(p) . Many inves- 
tigation have been done on the intermittent properties of 
this model @ and of the original GOY model 0|J,|o| 
at varying the coefficients such as to have different con- 
served quantities, in order to investigate the importance 
of inviscid invariants in determining the energy transfer 
properties. Some evidences supporting non-trivial effects 
introduced by non-positive invariants have been quoted. 
In this letter we discuss the dependency on the other im- 
portant free parameter entering shell-modelization: the 
separation between neighboring shells A. By decreas- 
ing A one describes turbulent energy transfer in terms 
of interactions which become more and more local in 
Fourier space. In the limit of A — > 1 one recovers a 1- 
dimensional partial differential equation j2^Jl^]. Local- 
ity of interactions has always been a long-debated issues 
of the K41 picture. In fig. 1 we show the 6th order 
structure function of model (||) for different inter-shell 
separations: A = 2 and A = 1.05. In both cases the total 
number of shells is chosen such that the physical length 
of the inertial range stays constant. Clearly, there is a net 
trend toward a less intermittent state by decreasing the 
shell-ratio. Nevertheless, the exceptionally good scaling 
properties allow us to estimate, using the Extended Self 
Similarity method ||, small deviations from K41 scaling 
also for the less intermittent value A = 1.05. 

In view of the discussion about log-poisson statistics 
this analysis could not be sufficient. Non- universal as- 
pects can be masked by trivial properties of the most 
singular structure. We have therefore tried to highlight 
possible universal aspects of the chaotic energy transfer 
by focalizing our attention on the non-linear part of in- 
termittent corrections and in particular on the parameter 
P which can be extracted by assuming log-poisson inter- 
mittency. The best way of extracting (3 from numerical 
data is to look if structure functions verify the log-poisson 
hierarchy 0: 



Sp+i(n) = A S p (n) > 
S p (n) " S p -i(n) 



(4) 



where f3' = /3 1 / 3 and A n has a scaling dependency on k n 
but it is p-independent. If structure functions follow a 
log-poisson statistics one can extract /?' from a linear fit 
of a log- log plot of (|i|), at varying p and fixed n. 
In fig. 2 we show relation (Q) for the two inter-shell sep- 
arations A = 2, 1.05 and for p = 1, 7. As it is possible 
to see, the straight line behaviour in the log-log plot is 
nicely verified (supporting the log-poisson assumption) 
and the slope of the two lines is the same, indicating 
that although the full intermittent corrections are differ- 
ent, the part directly affected by the (3 parameter is the 
same. In table 1 we collect all our estimates of f3' as a 
function of A. From these data one can safely conclude 
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that [3 stays constant by approaching the continuous limit 
A — ► 1. The remarkable change in the statistics reflected 
by changes in the ((p) exponents of fig. 1 may be only 
due to changes of the most singular energy-burst statis- 
tics, i.e., to changes of the parameters ho and do in the 
log-poisson jargon. Let us notice that already other au- 
thors have focused the attention on burst-like solutions 
in GOY models as possible explanation of its intermit- 
tent properties pjL53 |. Moreover, it is very simple to 
interpret the measured tendency toward K41 statistics 
as the results of a smoothing in the energy transfer due 
to the increasing local character of dynamical interac- 
tions, in the limit A — ► 1. In other words, we interpret 
the corrections to K41 for A >> 1 as the consequence of 
burst-structures travelling along the inertial range. Very 
singular bursts (with ho less than the Kolmogorov value 
1/3) appear only when there is a relevant mismatch be- 
tween eddy-turnover times of neighboring shells, i.e. only 
when shell-ratios are much larger than one. Otherwise, 
energy tends to flow smoothly toward small scales follow- 
ing K41 statistics. 

These bursts are non-universal in the sense that their 
statistics and their degree of singularity is a function of 
the inter-shell ratio, A. 

On the other hand, the mechanism underlying burst 
transfer seems only determined by the non-linear struc- 
ture of shell-model dynamics and it seems to fix a uni- 
versal value of (3 in the log-poisson description. 

To conclude, we have analysed and discussed the inter- 
mittent properties of a class of shell models for turbulent 
energy transfer at varying the inter-shell ratio parame- 
ter. We have shown that there exists a suitable limit, 
A — ► 1, for which intermittent fluctuations disappear and 
the scaling properties of the system become close to the 
K41 theory. By employing a log-poisson description of 
intermittent corrections, we have disentangled the limit 
A — > 1 in two main components: from one hand the most 
singular structure seems to tend continuously toward a 
K41 scaling (ho = 1/3); from the other hand the proba- 
bility distribution of intermittent fluctuations is still log- 
Poisson and described by the same defect parameter j3. 
One possible way of describing this result, in terms of 
the scaling exponents £(p), consists in saying that the 
quantity 



Ppq — 



C(p)- P /3C(3) 
C(9)-?/3C(3) 



(5) 



is constant for A — ► 1. Similar universal properties 
have already been reported in shell models at varying the 
dissipation mechanism, i.e. using hyperviscosities ]23j . 
Our results are qualitatively similar to what observed 
in real turbulent flows going from the inertial subrange 
to the viscous subrange. Indeed in this case anomalous 
scaling disappears (i.e. C(p)/C(3) —*p/3) while p pq stays 
constant, as recently observed in nw. 
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FIGURE CAPTIONS 

• FIGURE 1: Log- log plot of the 6th order struc- 
ture function vs k, for the two cases A = 2 (cir- 
cles) and A = 1.05 (diamonds). A linear fit in 
the inertial range, using Extended Self Similarity, 
gives: C(6)/C(3) = 1.983 ± 0.005 for A = 1.05 and 
£(6) /C(3) = 1.76 ± 0.01 for A = 2. 

• FIGURE 2: Log-poisson hierarchy (eq. ^) in a log- 
log plot for A = 2 (circles) and A = 1.05 (dia- 
monds). For each A, we took p = 1, 7 for three 
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different scales n in the inertial range and we shifted 
the sets along y-axis in order to perform a single 
linear fit (solid lines). 

TABLE CAPTIONS 

• TABLE 1: j3' at varying A. Each value is the aver- 
age of slopes evaluated as in fig. 2 for all scales n 
in the inertial range. The theoretical prediction |t]] 
is /?' = (2/3) 1 / 3 = 0.87. 



A 





1.05 
1.2 
1.5 
2 


0.86 ±0.01 
0.88 ±0.01 
0.88 ±0.01 
0.86 ±0.01 



